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Abstract 

We present the statistical-mechanical theory of semiflexible polymers based on the connection 
between the Kratky-Porod model and the quantum rigid rotator in an external homogeneous field, 
and treatment of the latter using the quantum mechanical propagator method. The expressions and 
relations existing for flexible polymers can be generalized to semiflexible ones, if one replaces the 
Fourier-Laplace transform of the end-to-end polymer distance, l/(fc 2 /3 +p), through the matrix 
P(k,p) = (I + ikDM)^ 1 D, where D and M are related to the spectrum of the quantum rigid 
rotator, and considers an appropriate matrix element of the expression under consideration. The 
present work provides also the framework to study polymers in external fields, and problems 
including the tangents of semiflexible polymers. We study the structure factor of the polymer, 
the transversal fluctuations of a free end of the polymer with fixed tangent of another end, and 
the localization of a semiflexible polymer onto an interface. We obtain the partition function of 
a semiflexible polymer in half space with Dirichlet boundary condition in terms of the end-to- 
end distribution function of the free semiflexible polymer, study the behaviour of a semiflexible 
polymer in the vicinity of a surface, and adsorption onto a surface. 

PACS numbers: 36.20.-r, 61.41. +e, 82.35.Gh 
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I. INTRODUCTION 



Polymers with contour length L much larger than the persistence length l p , which is 
the correlation length for the tangent-tangent correlation function along the polymer and 
is quantitative measure of the polymer stiffness, are flexible and are described by using the 
tools of quantum mechanics and quantum field theory . If the chain length decreases, 

the chain stiffness becomes an important factor. Many polymer molecules have internal 
stiffness and cannot be modelled by the model of flexible polymers developed by Edwards 

0. 

The standard coarse-graining model of a wormlike or a semiflexible polymer was proposed 
by Kratky and Porod p. A few first moments of G(r,N) were computed in The 
iterature on the earlier work on semiflexible polymers can be found in the book by Yamakawa 
11 1 . For recent work see (and the references therein). Despite a considerable interest 



and immense efforts in last decades there is no theory of semiflexible polymers providing a 
general tool for treating problems including semiflexible polymers. 

In this article we present the theory of semiflexible polymers based on the relation between 
the Kratky-Porod model and the quantum rigid rotator in an external field |4|, an d 
treatment of the latter in the framework of the quantum mechanical propagator method j3| . 
Although the most polymer quantities are defined through the positions of the monomers 
r ( s )? (0 < s < N) (the end-to-end distribution function, the scattering function, the isotropic 
monomer- monomer interactions, etc.), and on the contrary the quantum rigid rotator is 
formulated in terms of tangents t(s) = dr(s)/ds, we have shown that the relations for 
flexible polymers (polymers in external fields, polymers with self-interactions, etc.) can 
be generalized to semiflexible polymer, replacing the Fourier-Laplace transform of the end- 
to-end distribution function of the flexible polymer, l/(k 2 /3 + p), by the infinite order 
matrix P(k,p) = (I + ikDM)~ l D with matrices D and M related to the spectrum of 
the rigid rotator, and considering an appropriate matrix element of the matrix expression. 
The quantity P(k,p) plays the key role in the theory similar to the bare propagator in 
common quantum field theories. The end-to-end distribution function is simply the matrix 
element (o,0\P(k,p)\0,0^, the scattering function of the polymer is the inverse Laplace 
transform of G(k,p)/p 2 multiplied by 2/N, the partition function of the stretched polymer 
is Z(f,N) = G(k = —iF/ksT^N) etc. The elimination of summations over the magnetic 
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quantum number in intermediate states enables one to carry out the calculations of above 
quantities with the infinite order square matrix P s (k,p). The present theory provides also 
the framework to study semiflexible polymers in external fields and with self-interactions, 
and problems including the tangents of polymer configurations. 

The article is organized as follows. SecO introduces to the Green's function formal- 
ism of the quantum mechanical rigid rotator. Sec lIIII derives the exact expression for the 
Fourier-Laplace transform of the end-to-end distribution function, and establishes its im- 
portant properties. Sec JIVI presents results of the computation of the scattering function of 
a semiflexible polymer. SecfVl considers the localization of semiflexible polymers in a weak 
symmetric potential corresponding to adsorption onto an interface. Sec JVll considers the 
behaviour of the polymer in the presence of a surface. Sec JVllI treats the adsorption onto a 
surface. Sec JVIIll introduces to the description of semiflexible polymers with self-interactions. 



II. FORMALISM 



The Fourier transform of the distribution function of the end-to-end polymer distance of 
the Kratky-Porod model [f| G(k,L) = J d 3 r exp(— ik(r — r ))G(r — r ,L) is expressed by 
the path integral as follows 



G(k,L) = J Dt{s) \\5{t(sf~ l)expHk^ dst(s) - | J ds( 



(i) 



where l p is the persistence length, and t(s) = dr(s)/ds is the tangent vector at the point 
s along the contour length of the polymer. The 2nd term in the exponential is associ- 
ated with the bending energy. The product over s in Eq.Q takes into account that the 
polymer chain is locally inextensible. For a polymer which interacts with an external po- 
tential and for polymer with monomer-monomer interactions the terms — dsV(r(s)) and 
— | L ds2 Jq ds\U(r(s2) — r(si)) should be added in the exponential of Eq.((H), respectively. 

The path integral (0) (without the term depending on k) corresponds to the diffusion of a 
particle on unit sphere, |t(s)| = 1 J], lj . and is also equivalent to the Euclidean quantum 
25]. The Green's function of the quantum rigid rotator or that for diffusion 
of a particle on unit sphere obeys the following equation 

-^P (9, cp, L; O , ipo, 0) - ^V^P = 5(L)5(Q - Q ), (2) 



rigid rotator 



mi 
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where Q is the spherical angle characterized by angles 9, and <p, and 5(Q — Qq) is a two 
dimensional delta function having the property J dQ5(Q — Qq) = 1. Henceforth, instead 
of the contour length L we will use the number of segments N = L/l p . All lengths will 
be measured in units of the persistence length l p . The quantity P(9, (p, N; 9q, (po, 0) is the 
Fourier transform of the end-to-end polymer distance with fixed tangents of both ends. The 
Fourier transform of the distribution function of the end-to-end polymer distance is obtained 
from P(Q, N; Qq, 0) by integrating the latter over Q and Qq 

G(k, N) = ^- J dn J dtt P(tt, N; fi„, 0). (3) 

Notice that in the quantum mechanical counterpart of the problem iV corresponds to the 
imaginary time it. The bare Green's function Pq(9, (p, N; 9q, <po, 0) associated with (J2J) reads 

P (9, <p, N; Bo, <Po, 0) = ]T exp(- l{l + ^ )N )Y lm {9 ) V )Y* n {9 Q) <p ), (4) 

l,m 

where Yi m (9, <p) are the spherical harmonics, and I and m are the quantum numbers of the 
angular momentum. For a given I, m takes the values — /, — / + 1, I. The quantity 
Po{9, <p, N; 9q, ipo, 0) corresponds to Eq.Q with k = and with the following boundary 
conditions for the path t(s) (0 < s < N): t(N) = (9,<p), and t(0) = (9 ,<p ). 

We now will consider the Green's function P(9,ip,N;9o,<po,0) associated with Eq.fJTJ. 
The differential equation for P is 

-^P(9, <p, N; 9 , 0) - i V^P + U{tt)P = 6{N)6{Q - O ), (5) 

where CZ'(ktn) = ikt^ is the potential energy of the rigid rotator in an external field ik, 
where k is measured in units of As it is well-known from Quantum Mechanics j3| the 
differential equation (JHJ) can be rewritten as an integral equation as follows 

P(Q, N; fi„, 0) = P (Q, N] O , 0) - J ds j dQ'P (n, N; Q', s)U(ktcy)P(tf, s; Q , 0). (6) 

As we already mentioned above Eqs.(J3JE|) describes also the Euclidean rigid quantum ro- 
tator in an external field. The iteration of Eq.© generates the perturbation expansion of 
P(Q, N; Qq, 0) in powers of the potential [/(kt^), and can be symbolically written as 

P = P - PoUP + P UP UP - ... = P - P UP (7) 
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The coefficient in front of (k 2 ) n of the expansion of G(k, p) and consequently of P in powers of 
k 2 , multiplied by the factor (— l) r T(2n + 2) is the Laplace transform of the moment < r 2n > 
of the end-to-end distribution function. Thus, Eq. J7J) gives the moment expansion of the end- 
to-end distribution function. The integral equation (JBJ is nothing but the Dyson equation. 
The description of semiflexible polymers based on (JBJEJ) is a variant of the application of the 
methods of quantum field theory to problems of statistical mechanics j^j], ^| (and citations 
therein). 



III. THE END-TO-END DISTRIBUTION FUNCTION AND ITS MOMENTS 
A. General consideration 

We now will consider the computation of the end-to-end distribution function of a semi- 
flexible polymer using the propagator method. The matrix elements of the external potential 
U(Q) over the spherical functions are given by 

(l',m'\U\l,m)=ik [ dnY; m ,(6,<p)(nt n )YUe^) = ikM.^UOi^i), (8) 



where the unit vector t is characterized by the angles 9 and (p, and n = h/k by the angles Q\ 
and (f\. Due to the convolution character of the expression (JHJ) with respect to the integra- 
tion over the contour length (P Q (Q, N; Q f , s) depends on the difference N — s), the Laplace 
transform of P(fl, N; Qq, 0) in Eq.© with respect to N permits to get rid of integrations 
over the contour length. Thus, in the following we will consider the Laplace transform of 
G(k,N) with respect to N. Using the spectral expansion of Pq according to Eq.pjl in the 
perturbation expansion of P(Q, N;Qq,0) which is given by Eq.((7j) enables one to sum the 
moment expansion of G(k,p) in powers of k as 

G(M = (0,0|P(k,p)|0,0), (9) 

where 

P(k, p) = (I + ikDMY x D, (10) 
and D is infinite order matrix defined by 

A /' = T77-, 1 ; s i v (11) 
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with I, I' = 0, 1, .... Eq. ffHH) can be also derived via direct solution of the Dyson equation 
(JZJ). The zeros at the brackets in Eq.(J£JJ) means I = and m = 0. Notice that Eqs.(P?WTUJ) 
have to be understood as quantum mechanical expectation values, so that summations over 
the quantum numbers I and m of the angular momentum in the intermediate states of the 
series of P(k,p) are implied. 

We now will show that summations over the magnetic quantum number m in interme- 
diate states in Eq.Q can be eliminated. This can be shown using the fact that G(k,p) is 
function of k 2 , which legitimates us to choose the z-axes of the coordinate system along k in 
computing the matrix elements (/', m'\U\l, m). In this case the scalar product nt in Eq.(jHJ) 
becomes simply cos#, with the consequence that the matrix elements M^ >m / ; / jm become diag- 
onal with respect to indices m and m', so that the magnetic number will be zero throughout 
the products of matrices in the series of P(k,p). As a result Eqs.(p irU|) simplify to 



G(k,p) = (0\P s (k,p)\0) (12) 

with 

P s (k,p) = (/ + ikDM s Y x D, (13) 
where the square matrix M s is defined by 

Mf v = wtfiy+x + wi +1 6i +lj i,, (14) 

and wi = a// 2 / (4/ 2 — 1). Summations over the intermediate states in Eqs. (jl2lll3*j) occur 
over the eigenvalues of the angular momentum / = 0, 1, so that according to Eqs. (jl2HT??j) 
the calculation of G(k, v) reduces to the computation of the matrix element of an infinite 
order square matrix 23j] . Eqs.lJU HTUJ) and (|12H13|) . which are more general as those derived 
in 2s| , enable one to compute the end-to-end distribution function with fixed tangents, and 
to study polymers in external fields and with self-interactions. The validity of Eqs. ()12HT3~|) 
can be proved in more general way using the relation 

Zi,m 1 |P(k,p)|Z a ,m a ) =5>S mi (a,/3,7) (h\P s (k,p)\l 2 ) D^^a, (3, 7 ), (15) 



where D l m m ,{oc, /3, 7) is the Wigner D-function |2J, D l m m ,{a, (3, 7) = e- ima d l m rn ,e- irri '^ and 
a, f3, 7 are Euler angles chosen such that the z-axes of the transformed coordinate system 
is directed along the wave vector k. In obtaining (|15|) we have taken into account that the 
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matrix elements (li, m' x \P{k.,p)\l2, 77i 2 y computed in the coordinate system with the z-axes 
parallel to k are diagonal with respect to indices m' 1 and m' 2 . Using (JT5j) and the property 
of the Wigner D-function, D l mQ (a, f3, 7) = y 'An / (21 + T)Y* m {(3, a), one obtains 

(z,m|P(k,p)|0,0) = V4tt/(2^ + l)^i^i) (l\P s (k,p)\o) . (16) 

Eqs. ()15M16|l are valid as well for matrix elements of (DM) n , which are terms in the series 
of P(k,p) in powers of k 2 . Notice that summations over intermediate states on the left side 
occur over I and m, while summations on the right side are only over I, as it is already clear 
from notations. To establish Eqs. (|12llT5j) using we consider the expression 

(0,0|(L>M) ni+n2 |0,0) = ^2 (0,0\(DM) ni \l,m) (l,m\(DM) n ' 2 \0,0) . (17) 

The application of (|16j) to both off-diagonal matrix elements in ()17j) combined with 
the use of the addition theorem for spherical functions, Pz(cosa>) = 4tt/(21 + 
1) Ylm=-i Yim(0', y^O^mX^ V 9 ); with to being the angle between the vectors characterized by 
spherical angles 9', tp' and 9, ip, which is zero in the case under consideration, shows that the 
summation over m in (|17|) will be eliminated, and we obtain again Eq.(|12j). Eq.(|15|) enables 
one to replace P(k,p) in favor of P s (k,p) in the expression 

J2(h,m 1 \P(k 1 ,p)\l,m^l,m\P(k 2 ,p)\l 2 ,m 2 ^ = ^ D l £ >mi { ai , ft, 71) (k\P s (h,p)\l) 

X D mi,m( a l» A, 7l)^, m (a 2 , £2, 72) 

rn 

x (l\P s {k 2 ,p)\h) /J^ m2 (a 2 ,/5 2 , 72 )(18) 

where the Euler angles at, ft, 7 j define the coordinate system with the z-axes parallel to the 
wave vector kj. The limits in summations over the magnet quantum numbers in (|18jl are 
determined by the corresponding quantum number of the angular momentum. Eq. (jl8)l can 
be generalized in a straightforward way for matrix elements of products of arbitrary number 
of propagators P(kj,p). The sum over m in ()18|) can be carried out using the addition 
formula of the Wigner P-function 

J2 Dl m,m"( a h /5l,7l)^!n»,m'(«2,/5 2 ,72) = D l mm ,(a, 0, 7), 
m" 

where a, j3, 7 are the Euler angles for the resulting coordinate transformation S — » Si — >• S 2 , 
to give 

D l m>ml (a',(3,y'). (19) 
7 



The angles a 1 and 7' correspond to the coordinate transformation S — > S\ — > S' 2 with 
a' 2 = —02, = P2, I2 = ~l2- Expressions similar to (fTHj) appear in studies of semiflexible 
polymers in external fields (adsorption) and with self-interactions, where in the perturba- 
tion expansions of the quantities such as partition function, "propagators" P(kj,p) with 
different wave vectors k; appear. The consideration of the expression in the middle line of 
(I18|) containing the sum over m in one dimension, where the Euler angles take the values 
on = 7i = 0, f3i = 0, 7r which means that in this case the wave vectors k, can be only 
parallel or antiparallel. In this case the expression (fT§j) reduces to D l mm , (0,0,0) = b m ,m' or 
Dj„ m ,(0,f,0) = (— l) l+m 5 m - m f. Therefore, the magnet quantum numbers in the intermedi- 
ate states of expectation values of expressions like (|18|) over the ground state will be zero, 
so that the factors in the intermediate states (|19|) become simply (±1)'. The sign minus 
applies, if the neighbor wave vectors are antiparallel. 

In explicit computations based on Eqs. (|12llf3j) one should truncate the infinite square 
matrix P s (k,p) by a finite matrix of order n. The expression for G(k,p) obtained in this 
way is a rational function being an infinite series in powers of k 2 , i.e. it contains all moments 
of the end-to-end distribution function, and describes the first 2n — 2 moments exactly. In 
context of eigenstates of the rigid rotator, the truncation at order n takes into account the 
eigenstates with quantum number up to the value I — n — 1. The truncation of P s (k,p) by 
4 order matrix, which is the consequence of corresponding truncation of the matrices D and 
M s , yields the Fourier-Laplace transform of the end-to-end distribution function as follows 



1 1 + i£ + 

1 L ^ 15f14.nl (Z+r,) ' 



4fc 2 1 9fc 2 

15(l+p)(3+p) ~T 35(3 +p )(6+p) 



^p) = -—^- ;;r> n " " m — (20) 

y 1 ~T 3p(l+p) ~r 15 (l+p) (3+p) 35(3+p)(6+p) 35p (1+p) (3+p) (6+p) 

Ge{k,p) is the infinite series in powers of k 2 , and describes exactly the first 6 moments of 
the end-to-end-distribution function. The moments of the end-to-end distance are obtained 
from (fT2*|) as 

(R 2n - 2 ) = T(2n)C; 1 (0|GDM s ) 2n - 2 |0) , (21) 

where C~ l denotes the inverse Laplace transform with respect to p, which is the Laplace 
conjugate to N. Carrying out the inverse Laplace transform of (l2~Tj) using Maple or Mathe- 
matica we have analytically computed 28 moments of G(r,N) 23(. From the expansion of 
the exact formula (|12|) for large p at different orders of truncations we have obtained the 
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leading terms of G(k,p) as 



k ~\~ ~ k ~ k ~\~ ...^ 

p 3p 6 5p b 7p 7 



which is nothing but the series of the Fourier-Laplace transform G ro d{k,p) = k^ 1 arctan(A;/p) 
of the end-to-end distribution function of a stiff rod G ro d{r, N) = (4irN 2 )~ 1 5(r — N). Taking 
into account the subleading terms in the expansion of G(k,p) 

3p* k2 ~ 3^ + V &k& ~ W* k% + "' 
results in the following expansion of the end-end-distribution function for small N 

N d 

G(r, N) = G rod (r, N) - -—G rod (r, N) + ... (22) 

aJM 

The latter can be considered as a singular expansion of the end-to-end distribution function 
over its width. In principle, the expansion ()22|) can be extended to take into account the 
next terms. However, so far we could not sum the expansion of the next to subleading terms 
in Eq.(fT2i) 

4(4(-) 2 + t|(-) 4 - 1(-) 6 + f A 8 - t(-) 10 + ^r(-) 12 - -)• 

p 6 6 p 15 p 15 p dp dp dp 
B. The spectral representation of G(k,p) 

The spectral expansion of the matrix A = DM enables one to derive the representation 
of P s (k,p) through the eigenvalues and the eigenvectors of the matrix A. The latter is 
useful in the treatment of polymer adsorption, which is considered in SectJV] The matrix A 
is not symmetric, so that to find its spectral expansion one should consider the eigenvalue 
problem for both A and the transposed matrix A. The eigenvalues of A and A coincide, while 
their (normalized) eigenvectors are different and are denoted by u(n) and u(n), respectively. 
Thus, the spectral decomposition of A reads 

Aij = 2J \ n Ui{n)uj(n) . (23) 

n 

Using (J2~3~j) we obtain P s (k,p) as 

[l\P s (k,p)\l') = 1+ ) kX m(n)Tfr(w)Ay'- (24) 

n n 



The eigenvalues X n of A truncated by even order matrix build pairs with the same absolute 
values and opposite sign in the pair, which we denote by ±Afc. For odd order A one eigenvalue 
is zero, and the remainder build pairs similar to those for even order matrices. The largest 
eigenvalue of A approaches the value l/y/Sp for p — > (the flexible limit). These properties 
of the eigenvalues of A enable one to write G(k,p) as 

G(k,p) = (0\nk,p)\0) = Lj2 TT ^l m , (25) 

where we have introduced the notation l m = u (m)u (m). The summations in (}25|) has to 

be carried out over all eigenvalues. We explicitly computed the eigenvalues and eigenvectors 

of A with Maple or Mathematica using truncations until n — 8. The eigenvalues and the 

factors l m using truncation with four order matrices are obtained as 

- _ , — (p(18 + 27p+ 10p 2 +p 3 )(85p+ 15p 2 + 105 + A)) 1/2 
1 ~ 630p + 945p 2 + 350p 3 + 35p 4 

_ (p(18 + 27p + 10p 2 + p 3 )(85p + I5p 2 + 105 - A)) 1 / 2 
3 ~ 630p + 945p 2 + 350p 3 + 35p 4 

with A = (7540p 2 + 1500p 3 + 15960p + 120p 4 + 11025) 1 / 2 . 

h = ( 3A " 10 p 2 + 6 °p + 315 )' z 3 = tJa ( 3A + 10 ^ 2 " Q °P ~ 315 ) 

The use of the spectral representation of G(k,p) given by Eq.()25|) enables one to carry out 
easily the inverse Fourier transform of G(k,p), and reduces the computation of G(r,N) to 
performing the inverse Laplace transform of G(r,p). However, the computation of G(r,N) 
using truncated expressions results in G(r, N) taking negative values for large r, and de- 
mands a special consideration. Very recently the computation of G(r, N) by using different 
methods was considered in [3|-[l^. 

Notice the following difference in truncation of (piTUjl with odd and even order matrices. 
While G(k,p) for even n behaves for large k as 1/k 2 , it behaves as const using truncations 
with odd n. The moment expansion as well the expansion for large p of G(k,p) behaves 
correctly. 



C. The Markovian property of the end-to-end distribution function 



It is well-known that the end-to-end distribution function of an ideal continuous flexible 
polymer (and of a polymer in an external potential) possesses the Markovian property, which 
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takes in terms of the Fourier transform Po(k, N) = exp(—k 2 N/d) of the end-to-end distance, 
the simple form, Po(k,N) = Po(k,N — s)Po(k,s). The Markovian property is obviously 
not valid in this form for a semiflexible polymer. Nevertheless, it can be generalized in 
an appropriate form to a seminexible polymer, too. Using the definition of the end-to-end 
distribution function (|1J) we obtain in a straightforward way 



G(t-t ,N) = J ^exp(zk(r - r ) - ik dst(s)) 

= J d s r' /j exp(zk 2 (r - r') - ik 2 jf dst(s)) 

x f exp(zk x (r' - r ) - iki f dst(s))\ , (26) 

Jkl JO / 

where s' fulfils the condition < s' < N . Representing the average in (J26|) over the 
eigenstates of the rigid rotator we arrive at 

(0\P s (k, N)\0) = (Q\P S (k, N - s)\l) (l\P s (k, s)\0) , (27) 
i 

where P s (k, N) is the inverse Laplace transform of the matrix P s (k,p). Eq. (j27|) is the gener- 
alization of the Markovian property for a seminexible polymer. Notice that the summation 
in (|2*7jl occurs only over the quantum number of the square of the angular momentum /. Sim- 
ilarly, one can show that the Markovian property applies to the off-diagonal matrix elements 
(l\P s (k, N)\l'), too. Eq . (|27p can be immediately generalized by partitioning the interval (0, 
N) in Ss = N/n intervals. The latter is an important ingredient of our consideration of the 
behaviour of a seminexible polymer in half space (see Sec lVI|) . 

Notice that the Markovian property of the quantity (0, 0|P(k, iV)|0, 0) 

(0, 0|P(k, iV)|0, 0) = (°> °l P ( k ' N ~ s ) l Z ' m ) & m \ P ( k ' s ) 1°' °> ( 28 ) 

is the consequence of the group property of the time evolution operator of the rigid rotator 
in an external field. Applying (fTBf) to both transition amplitudes in (J2EJ), and using the 
addition theorem for spherical functions gives again Eq. (J27|) . 
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D. Tangent correlations 



We now will consider two examples of computation of quantities containing tangents. 
The correlation function of tangents (t(s)t(s')) can be calculated using (J3J as 

(t(s)t(s')) = Yl Yl < ' °l U ^ m > e-^(' +1 )l s - s 'l (I, m\ t t |0, 0) = e^" 5 ' 1 , 

i=x,y,z l,m 

where the arc length is measured in units of l p . 

The use of Eq. fll6|) enables one to derive the following exact expression for the Fourier 
transform of the end-to-end distribution function with the fixed tangent t(0) 



G(t,k;iV) = 5^(0 | P s (k,N) | /) y/^±ip,(tn). (29) 

Eq . (|29|) enables one to compute the transversal moments (with respect to the direction t(0)) 
of the free end of the polymer r(N). In the case if the tangent t(0) is parallel to the z-axes, 
the scalar product tn in Pj(tn) in (}2*9*|) becomes zero. The calculation of the 2nd and the 
4th transversal moments using Eq.(|29p in the case if t(0) |T e z yield 

<^> = ^(6iV-e-^ + 9e^-8), 

(l A v _ 8 2 _ 248iV 1834 J_ _ J_ 

\ tr) 9 l 15 + 75 + 1225 25 + 

S0_^ N _624 e _ N + ±Ne _ 3N _3 ANe _ N) 



147 25 21 5 

The computation of higher moments is similar. For small N the moments behave as (rf r ) ~ 
(2/3)iV 3 , (rf r ) ~ (8/9)iV 6 . The ^-dependence of (rf r ) on N means that the transversal 
fluctuations, which are controlled by the bending energy, are small in comparison to the 
length of the polymer N. For large N the moments behave as (rf r ) ~ (4/3)iV, (rf r ) ~ 
(32/9) A^ 2 , and obey the relation between the moments of a Gaussian distribution in two 
dimensions, (rf r ) = 2 (r 2 r ) 2 . Note that his relation is also fulfilled for small N. The quotient 
i r tr) I i r tr} 2 nas the minimum equal to 1.36 at iV = 1.82. 

IV. THE STRUCTURE FACTOR OF A SEMIFLEXIBLE POLYMER 

In studying the structure factor of a semiflexible polymer which is defined by 

2 pN pS2 

S(q, N ) = ^J o ds * J o ds ± (exp(zq(r( S2 ) - r( Sl )))> . (30) 
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we first express r(s 2 ) — r(si) in (|3*0|) through the tangent vectors, r(s 2 ) — r(si) = dst(s). 
Representing the average in (j3Tlj) through the eigenstates of quantum rigid rotator yields that 
the structure factor of the semiflexible polymer S(q, N) is the inverse Laplace transform of 
G(q,p)/p 2 multiplied with the factor 2/N J^. Fig^ shows the double logarithmic plot 
of the structure factor of a semiflexible polymer as a function of the absolute value of the 
scattering vector q (measured in units of l p ) using truncations of the exact matrix expression 
with finite order matrices. The slope —1, which is characteristic for rigid rod behaviour, is 



N=0.7 



w 

CO 

o 



Slope -1 




log[qi p ] 



FIG. 1: The log- lot plot of the structure factor of a semiflexible polymer. The dashed curves: 
truncations with number of exact moments n = 14; continuous curves: n = 18. 



also shown as a guide for eyes. Note that the curves for n = 8 and n — 10 takes exactly into 
account 14 and 18 moments of the end-to-end distribution function, respectively. 

The structure factor of a semiflexible polymer approaches for small N that of a stiff rod. 
The latter can be computed using G(k, N) = sin(giV) / qN , carrying out the inverse Laplace 
transform of G(k,p)/p 2 , and multiplying it with 2/N. The result is 

2 



S(q,N) rod 



Nq- 



: (NqSi(qN) + cos(qN) - 1), 



where Si (or) = J* dtsxn(t)/t. The subleading terms in the expansion of G(k,p) for large p 
given by (J22)) result in the following correction to the structure factor of the stiff rod for 
small iV 



2 sin(giV) cos (q N) 



+ 



q 3 N ' 3q 2 
The plots of qS(q, N) and qS(q, N) rod are shown in FigJ21 



(31) 
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FIG. 2: The structure factor of a semiflexible polymer, qS(q). Dashed curves: n = 14, continuous 
curves n = 18. The structure factor of rigid rod plotted for N = 0.7 coincides practically with the 
dashed curve. 

The computation of the structure factor is exact until the values of q where the curves 
associated with different truncations begin to diverge. Thus, the present method enables 
one in fact an exact computation of the structure factor of the Kratky-Porod model without 
restriction on the polymer length. The structure factor computed in (2^ describes the 
low q range corresponding to the Gaussian coil region. The exact structure factor for the 
infinitely long chain given in 



28| is valid only for large q. The approximative approach used 
in [31] does not go beyond the exact description of the exact 2nd moment of the end-to- 
end distribution function. The expression of the structure factor obtained in j^] gives an 
interpolation between the Gaussian coil and stiff rod limits. 

Besides the interest in its own as an experimentally accessible quantity, the structure 
factor can be used within the random phase approximation P to study the effects of rigidity 
on the phase behaviour of polymer mixtures of different architectures (3^]. 



V. ADSORPTION IN A WEAK SYMMETRIC POTENTIAL 

We now will consider the adsorption of a semiflexible polymer in an external delta- 
potential U(z) = —u5(z), which corresponds to adsorption onto an interface placed at 
z = 0. The strength of the potential u is measured in units of ksT /l p . First we consider the 
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partition function of the polymer with ends fixed at r and r' in a general potential U(r), 
which reads 

Z(r,N;r') = /$(r - r(JV))<y(i^ - r(0)) exp(- jf dsU{r{s)))\ . (32) 



It can be shown in a straightforward way that the Taylor series of (I32J) in powers of the 
interaction potential can be written as 

Z{t,N-y') = jT^exp(ik(r-r / )-ikjT tfet(s))^ + XJ(-l) n J d 3 r n ... J d 3 ri 



x ... C/(q n )...C/(qi) / ... / exp(iq n r n + ... + zqiri 
J q n J qi J fcn+i >^ fci 

-Hk n+1 (r - r n ) + ik n (r n - r n _i) + ... + iki(ri - r')) 
x( / ds„ / dsn-i— / <isiexp(-ik n+ i / dst(s)) 

\io JO Js n 

x exp(— ik n / rfst(s))... exp(— iki / dst(s)) ). (33) 

J s n -i Jo I 

Performing the integrations over r$ and kj establishes the equivalence of Eq. (j33p with the 
Taylor series of (J32j) . Expressing the average in (J33j) through the eigenstates of the rigid 
rotator, and carrying out the Laplace transform with respect to iV we obtain the bracket 
in the 2nd term of ()33|) as 



(...) = (|o,0|P(k n+ i,p)|/ n ,m n ^/ n ,m n |P(k n ,p)|/ 

n— 1) ^n— 1 

'z 1 ,mi|P(k 1 ,p)|0,0\, (34) 



where p is Laplace conjugate to N, and P(k n ,p) is given by Eq.(fTU|) . The summations over 
li and rrii (i = 1, n) are implied in ()34|1. 

The partition function Z(z,p; z') of the polymer in an adsorbing potential ?7(z) = — u8(z) 
can be obtained from Eq. ()33|) integrating it over r", so that k^ 1 in ()34j) become zero. The use 
of Eqs. (jl8|) permits to replace the matrices P(h n ,p) in favor of the square matrices P s (k n ,p). 
As it is shown in Sect 1IIIA1 the factors (±1) *, where the sign minus corresponds to the case 
when kf and k*_ l have different sign, appear in the intermediate states. It can be directly 
shown that the factors (— can be taken into account by allowing ki in P s (ki,p) to take 
positive and negative values. Inserting the delta-potential U(z) = —u5(z) into Eq.(j22I) we 
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obtain the partition function of adsorbed polymer as 

Z(z,p;z') = exp(ik{z- z'))(o\P s {k,p)\0^ 

xu (1x1(1 -(«/2tt) y dkP s (k,p))-% 

x jf expHfez') ^ 2 |P s (A; 2 ,p)|0) , (35) 

where I is a diagonal matrix of infinite order. The poles of the partition function are given 
by the zero points of the determinant 

/oo 
dkP s (k,p)) = 0. (36) 
-oo 

Inserting P s (k,p) = (I+ikDM s )~ l D = (I+k 2 (DM s ) 2 )- l (I-ikDM s )D into (JHED we obtain 
the eigenvalue condition in the form det(I — (u/2tt) dk(I + k 2 (DM s ) 2 )~ l D) = 0. We 
would like to stress that Eqs. (|35H3*Bj) are exact for adsorption in a Dirac's delta-potential. 
The formula for the flexible polymer is obtained from Eq. (j33j) by neglecting the off diagonal 
matrix elements and using the expression P s (k,p) = l/(k 2 /3 + p) for the propagator. 

It is interesting question if the energy eigenvalue condition (|36|) can be rewritten in 
terms of the boundary conditions imposed on "a wave function" as it is the case for a 
flexible polymer. Such an interpretation of Eqs . (|35ll3l)|) would enable one to consider the 
adsorption in potentials with finite widths and depths. An heuristic attempt to treat the 



22| . Notice that according to the analogy of 



adsorption in this way was undertaken in 
(|33|) with the corresponding equation for a flexible polymers it is tempting to interpret 
J ki exp(ikiz) (o\P s (ki,p)\lij as the "wave function" of the localized semiflexible polymer. 

The integrations over k in (|3~K|) can be carried out using the spectral representation 
according to 

/ exp(±tkz)^l 1 \P s (k,p)\l 2 ) =Y,fm\ z ) u h( m )uh(™>)D hth . 

with ft\z) = f,exp(±ikz){l - ik\ m )/(l + k 2 X 2 m ) = {Xm 9 l ± Xm exp(- \z/X m \). Note that the 
integral dzfm (z) is equal to 1. The inverse Laplace transform of ()35j1 . which takes into 
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account for large N the main contributions associated with poles defined by Eq. (j3Hj) . yields 

Z(z,N;z') bound ~ ^^^/^(^Hs^mjA,,!, 

po h,h m 

A h,i 2 (po) 



X ■ 



where 



^2f^(z')u l2 (n)uo(n)D , , (37) 



/oo 
dkp s (k,p))/(p-p ), 
-oo 



and Ai lji2 (p ) is the adjoint of the matrix element (1x1(1— (u/2tt) dkP s (k,p))~ 1 \l 2 
taken at po, which is the zero of (|36|) . Eq. (|37j) is the spectral expansion of the partition 
function over the localized states, and shows that the zeros po of Eq. (|3T)Jl yield the energy 
spectrum of the localized semiflexible polymer. The distribution function for monomers n(z) 
is calculated using the expression 

H ^ = fdzfdz^Z(z N-z>) J dS j dZ ' I dz " Z ( z '> N ~ s i z \i z ( z , s;z")i,o, (38) 

which generalizes the corresponding formula of flexible polymer Q]. The quantity 
Z(z', N; z)oj is obtained from the expression on the right-hand side of (|3*7)l by replacing 
the index in Uo(n) and in D 0tQ by I. Further, in computing n(z) we will take into account 
only the ground state (ground state dominance). Inserting (|3~7jl into (J3~%j) results in the 
following expression for the monomer density of adsorbed polymer 

n ( z ) = Yl Yl T7^ A v^{Po)fi~X z ) u h{ n )ui{n)D l rf£\z)u l ^ (39) 

n,m l 2 ,l,h ^°' ) 

The distribution function of one polymer end of adsorbed polymer, which is defined by 

f{z) = [ dz'Z(z,N;z')/ [ dz [ dz'Z(z, N; z'), 



results after inserting (|3*7j) and using the approximation of ground state dominance in the 
following expression 

f(z) = Yfn + \ z )Mn)ui(n)D l:l A lfi (p )/A 0fi (p )/D 0fi . (40) 

l,n 

We have mentioned in Sect llll S\ that G(k,p) behaves differently for large k in using 
truncations with odd or even order matrices. This is now important because Eq. (JBTIj) contains 
integrals over k of P s (k\,p). Notice that the matrix elements in (j3*Uj) which are odd in k 
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vanish because the integration is carried out in symmetric limits. While G(k,p), which 
corresponds to an odd truncation, does not have well-defined inverse Fourier transform with 
respect to k due to the divergence of the integral for large k, we will use truncations with 
even order matrices. At present we do not have a more reasonable explanation for this 
behaviour of G(k,p) for odd and even n for large k. 

Thus, we will use in Eq. (j36|) truncations with even order matrices, and will study the 
energy eigenvalue condition for n = 4 and n — 6, where 6 and 10 moments of the end-to-end 
distribution function are exact, respectively. The use of the spectral decomposition of the 
matrix DM S according to (|23j) permits to perform easily the integration over k. As a result 
we have obtained from (|36|) the following relations between the adsorption strength u and 
the localization energy E in the vicinity of the localization transition 

u = ^TE- l^(-E) - v / 3^("^) 3/2 + 0{{-E)% n = 4 (41) 

u = -^=V^E- 2.264(-£) - 2.0854(-£) 3/2 + 0((-£) 2 ), n = 6 (42) 
V3 

The 1st term on the right side of Eqs. (j41H4*2"j) is the result for a flexible polymer written in 
units under consideration, where the dimensionless adsorption strength and the energy are 
measured in units of ksT/lp and l p , respectively. The relation I = y/2l p allows to replace 
the persistence length in favor of the statistical segment length I. Eqs. (141114*2")) show that the 
applicability of adsorption theory of flexible polymer is in fact restricted to the localization 
transition u — > 0. For all finite u there are corrections due to the polymer stiffness. 

The normalized monomer density at the localization energy —E = p = 0.8, which is 
computed from Eq.(JHHJ) using truncation with matrices of order 4, reads 

n{z) = 0.42e" 3 - 842 + 2.35e- 12 ' 78 * + 4.52e" 21 - 732 . (43) 

In contrast, the monomer density for a flexible polymer for the same localization energy is 

n fl (z) = 1.55e- 3 - 10z . 

Eq. (J43|) shows that the decay of the monomer density of adsorbed semiflexible polymer is 
not exponential. This is also seen in FigOl displaying the logarithmic plot of the monomer 
density for two different localization energies. The monomer density for flexible polymer 
is also shown for comparison. As it follows from Eq. ()39|) and from the expression of the 
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FIG. 3: The normalized monomer density of adsorbed polymer. Solid line: adsorption energy 
p = 0.8, adsorption strength u = 0.375; dots: p = 0.3, u = 0.305; dashes: flexible polymer at 
p = 0.3, u = 2p/ v / 3 = 0.346. 

factor fm(z) the decay of n(z) away from the interface is controlled by the length |Ai| /2, 
where Ai is the largest eigenvalue of the matrix DM S . At z smaller than z c « 0.22, which 
is determined by the condition that the first term in (}4*3*j) have the same value as the rest, 
the decay of n(z) is sharper, and is determined by the length | A1A2I /(|Ai| + [ A2 1 ) , where A2 
is the next largest eigenvalue of DM S . The behaviour for z < z c is likely to be ascribed 
to formation of the liquid-crystalline ordering of the polymer caused by alignment of the 
pieces of adsorbed polymer along the interface. The length z c separating the faster and the 
slower decays, can be interpreted as the correlation length of the liquid-crystalline ordering 
induced by the interface. The analysis of the behaviour of n(z) for n = 4 shows that in 
approaching the localization transition po — > the prefactors in front of the 2nd and 3rd 
exponents in Eq . (|43|) tend to zero faster than that in front of the first term, which tends 
to the value y/3po, which is the result of the flexible polymer. For truncations with n x n 
matrices the monomer density n(z) is according to Eq. (j39|) a superposition of n(n + 2)/8 
exponents. We expect that similar to the case n = 4 the monomer density tends to that of 
the flexible polymer in approaching the localization transition, p — » 0. 

It is surprising that the leading correction in Eqs. (J41M2j) to adsorption of flexible polymer 
depends on the order of truncations. We expect that at given adsorption strength there are 
pieces of the polymer with characteristic size depending on u, which are completely adsorbed 
and are aligned along the interface. The order of truncation should be such that the statistics 
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FIG. 4: The normalized distribution of one polymer end of adsorbed polymer. Solid line: adsorption 
energy p = 0.8; dots: p = 0.3; dashes: flexible polymer at p = 0.3. 



of these pieces is described accurately. With increasing u the size of these pieces increases 
too, which demands more accurate description by using the higher order matrices. However, 
the extreme narrowness of the delta-potential demands high accuracy on quite small scales, 
which is apparently the reason why the leading correction is different for n = 4 and n = 6. 
Fig0] shows the distribution function of one polymer end of adsorbed polymer computed 
using Eq. ()40|) . Similar to FigJ3]we see here different dependence on z for small and large z. 

Note that truncation of Eq.(|36|) with 2x2 matrices, where only the 2nd moment of the 
end-to-end distribution function is taken exactly into account, yields the leading correction 
to the flexible limit as (— E) 3 ^ 2 instead of (— E) in accordance with Eq. (j4*T| . Thus, truncation 
with n = 2 does not correctly describe the leading correction to the flexible polymer. 

Eq. (j41)) shows that the correction term is negative, i.e. the localization energy for semi- 
flexible polymer at the same strength of the attraction potential is hig her than that for 
flexible one. This conclusion is in agreement with preceding studies jsj Adsorption 
of a semiflexible polymer considered in 40|-Q|, was carried out in the framework of the 
model, which is semiflexible on all scales. This model cannot describe the regime of weak 
adsorption considered in the present work. 
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VI. A SEMIFLEXIBLE POLYMER IN HALF SPACE 



The Fourier transform of the end-to-end distribution function (|T2*j) can be written in 
terms of the cumulants dsjj,2 n (s) of the moments as 

oo r N 



G(k,N) = exp(- V / d Sf i 2n (s)(k 2 ) n ) 

n=l J( > 



It is easy to see that the latter results in the following differential equation for the end-to-end 
distribution function 

dG{r dN° ,N) - B-l) n+ V 2 n(iV)A?G(r - ro, N) = 0, (44) 

n=l 

where A = V 2 is the Laplace operator. Eq.([44p. which is exact, enables one to derive the 
partition function of a semiflexible polymer in half space z > with Dirichlet condition at the 



boundary Taking into account that G(r = a/ (r" — i-q) 2 + {z + zo) 2 , N) also obeys Eq.(jHJ), 
the partition function of the semiflexible polymer in half space with Dirichlet boundary 
condition at z = can be written as 



Z(r, r , N) = G(^(r»-r<>) 2 + (z-z )\ N) - G(yJ(r» - r^) 2 + (z + z )\ N). (45) 

However, it is not clear from the above derivation if the trajectories of the polymer associated 
with Z(r, r , N) obey the Dirichlet boundary condition at intermediate points, < s < N. 
The following derivation of ()45|) based on the use of the Markovian property of the end-to- 
end distribution function enables one to prove this. Partitioning the interval (0, N) in n + 1 
equal intervals 8s = N/ (n + 1) we obtain 



G{t-t q ,N) = [ d\ n ... ! d 3 n ! ... / 



X 

exp(ik n+ i(r - r n ) + ik n (r n - r n _i) + ... + ikxfa - r )) 
(0\P s (k n+l , N - s n )P s (k n} s n - s n - 1 )...P s (k 1 , Sl )\0) . (46) 

Notice that ki in the argument of P s (ki,Si — Sj_i) is the absolute value of kj. The 
dependencies on positions r m in Eq. (j46|) are the same as in the corresponding equa- 
tion for a flexible polymer, which is obtained from ()46|) replacing the factors P s (k,p) by 
Pfi(k, N) = exp(—k 2 N/3), and taking into account in the intermediate states only the term 
with 1 = 0. In order to obtain from (J46)) the partition function of the polymer in half 
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space we will proceed in the same manner as for flexible polymers, and replace the factors 
exp(ik m (r m — r m _i)) in (j4T)j) for m = 1, n according to 

exp(ik m (r m - r m _i)) -> - exp(zk^(r^ - r^)) 

x (exp(i(k z m (z m - - exp(i(k^(z m + z m _a)))). (47) 

For m = n + 1 we use the replacement (|4"7j) without the factor 1/2. Inserting the latter 
into Eq. (|46j) gives the partition function of the semiflexible polymer with fixed ends in the 
presence of a surface at z = as 

Z(r,v ,N) = I d 3 r n ... I dh x I ... I exp(^ +1 (r ll -<)+zk!„(r;;-r^ 1 ) + ... 

+ik{(rl - r'^))2- n {exp{tk z n+1 {z - z n )) - exp{ik z n+1 {z + z n )))... x 
(exp(ifc*(zi - Zq)) - exp{ikl(z x + z ))) x 

(0|P s (A; n+1 , N - s n )P s (k n , s n - s n ^)...P s (k 1: Sl )|0) . (48) 

Z(r, tq, N) becomes zero, if z, z or any intermediate coordinate z m (m = 1, ...,n) are zero. 
Thus, the expression ()48|) for Z(r, r , iV) fulfils in the limit n — > oo the Dirichlet boundary 
condition z(s) = (0 < s < iV), and consequently Eq.(|48|) gives the partition function of 
the semiflexible polymer with both ends fixed in the presence of a wall at z — 0. 

We now will show that in the limit 5s — > Eq.()48|) passes over to Eq. (|45|) . To arrive at 
this result we carry out successively integrations over r m in ([48)1 . The integration over 
gives (27r) d ~ 1 5(k^ n — k^ +1 ) while the integration over z m yields 

27r(exp(i/c m _ ) _ 1 2 m _(_i ik m z m -i)5{k m ^ k m ) exp(ik m ^_ x z m +i + ik m z m -i)5(k m+ ^ k m ) 
exp(ik z m+l z m+1 - %k m z m _ x )b{k z m+x + k z m ) + exp(ik z m+l z m+1 + ik m z m _ x )b{k z m+x + k z J) (49) 

In obtaining (}4T?j) we have taken into account that P s (k m , s m — s m -i) depends on the absolute 
value of k m , and consequently can be written as common factor. Using the Markovian 
property of the distribution function (j2*7j) at the interval (s m +i, s m -i)> we obtain again the 
expression (|48|) with the difference that the point s m is now missed, and the prefactor in 
front of (}4*%j) will be 2~ n+1 . Repeating this procedure n times we finally obtain 

Z(r,r ,N) = I exp(ik ll (r" - r^)) {ex V {ik z {z - z )) - exp{ik z {z + z ))) (0\P(k, N)\Q) 
Jk 

= G(y^ - r )2 + ( Z - Zq )\ N) - G(yJ(r» - r'^ + (z + z ) 2 , N). (50) 
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The above derivation guarantees that the trajectories of the polymer r(s) contacting the 
surface at an arbitrary point s G (0, N) do not contribute to Z(r, tq, N). 

The derivation of (|48H5U|) is based on the Markovian property of the end-to-end distri- 
bution function given by Eq. (J28|) . and the dependence of P s (k, N) on the absolute value of 
the wave vector k. The statistical weights of configurations, which are determined by the 
expression in the brackets in Eq. (J4"%|) and are different for flexible and semiflexible polymers, 
do not play a direct role in the above derivation. The circumstance that the subtracted 
terms, which result in the 2nd term of Eq.([50jh obey the infinite order differential equa- 
tion of the free polymer (|44jl . guarantess that the polymer configurations contributing to 
Z(r, r , N) are those of the free polymer. Note that Eqs. (|45l5U|) are also valid for stiff rod, 
|r — r | = N, where the second term, which is not zero only if z or z are zero, selects the 
conformations, which do not have contact with the boundary. 



A. Behaviour of a semiflexible polymer in the vicinity of a wall 

Eq.(|50|) enables one to study the behaviour of a semiflexible polymer in the vicinity of a 
wall. For this goal we will use the approximative but simple analytic form of the end-to-end 
distribution function derived in [16] 

1 9L 
Go(r, L) = j^—^ e*p(- gai _ r2) ), (51) 

where r in Eq. (|51|) and throughout this subsection is measured in units of the contour length 
L, so that as a consequence that the polymer chain is locally inextensible, the distance r 
fulfils the inequality r < 1. The computations can also be carried out with end-to-end 
distribution functions derived recently in [3]>0I- 

To compute the distribution function G(r, r , N) with one end of the polymer fixed at 
distance z to the wall, we integrate Eq.(|50|) over r" — rJJ and zq, and obtain 

w(z) = 7r / = / dxg (x), (52) 



Jl-z 2 VI — P JO 

where go{x) is the normalized Gq(t,L) and x = 1 — r 2 . The integration over x can be 
performed analytically, while the integration over (3 only numerically. FigEl shows w(z) 
for different ratios L/l p . For z — 1 the polymer just begin to contact the wall, so that 
w(z = 1) = 1. For large L/l p the semiflexible polymer is a coil with the size which is 
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proportional to yL, so that the effect of the wall on w(z) appears at smaller z. For small 
L/l p the polymer behaves as stiff rod, and w(z) decreases linearly at small z. FigEl shows 
that at given distance of the polymer end to the wall z, and for given contour length L, the 
distribution function w(z) is larger for polymers which are more flexible. 
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z/L 



FIG. 5: The distribution function of one polymer end as a function of the distance to the wall for 
different persistence lengths. Solid line: L/l p = 10; dashes: L/l p = 1.5; dots: L/l p = 0.3. 

To compute the projection of the mean-square end-to-end polymer distance parallel to 
the wall under the condition that one end of the polymer is fixed at distance z to the wall 
we average (r" — Tq) 2 using Eqs. (j5()H5*TJ) . The normalization of the result with w(z) gives 

<(r"-ri5) 2 )^ = 7r f ~^=f , dx(J3-x)g (x)/w(z). (53) 

Jl-z* Vl- P J0 

FiglHl shows the mean-square end-to-end distance parallel to the wall as a function of z for 
different values of L/l p computed using Eq.(|53)l. While for Gaussian polymer the depen- 
dencies on longitudinal and transverse distances separate, ((r" — i"o) 2 ) z does not depend on 
z. In rod limit, L/l p <C 1, ((r" — ro) 2 ) 2 is a linear function having the values 1 and 2/3 
at z = and z = 1, respectively. The approximate distribution function (JSTj) reproduces 
qualitatively the rod behaviour. 

Fig E| shows the force f(z) = d\nw(z)/dz acting on the polymer as a function of the 
distance of one polymer end to the wall. As in the case of the distribution function w(z) 
the force begins to deviate from zero for more flexible polymer at smaller z. For very stiff 
polymer there are two regimes in the behaviour of the force. Just below z < 1 the force 
increase is sharp, and becomes weaker with further decrease of z. We attribute this initial 
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FIG. 6: The in-plane mean-square end-to-end distance of the polymer as a function of the distance 
of one polymer end to the wall for different persistence lengths. Solid line: L/l p = 10; dashes: 
L/l p = 1.5; dots: L/l p = 0.3. 

sharp increase of the force to the 
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FIG. 7: The force acting on the wall as a function of the distance to the wall for different persistence 
lengths. Solid line: L/l p = 10; dashes: L/l p = 1.5; dots: L/l p = 0.3. 

just begin to contact the surface, z < 1, the bending gives the major contribution to the 
force. With further decrease of z the configurations with large bending are less probable, 
and the increase of the force is of entropic origin. 



bending of the rod. It is intuitively clear that if the polymer 
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VII. ADSORPTION OF A SEMIFLEXIBLE POLYMER ONTO A SURFACE 



In considering the adsorption of a semiflexible polymer in a weak surface potential defined 
by U(z)) = —u5(z — Zq) for z > and U(z)) = oo for z < 0, where the attractive delta- 
potential is placed at the distance z to the wall, we chose as a reference state the distribution 
function in half space (|50p. and expand the partition function in Taylor series in powers of 
the attractive part of the potential. As a result we obtain the Laplace transform of the 
partition function of the polymer Z(z,p; z') as 

Z(z,p;z') = J {exp(ik{z - z')) - exp(ik{z + z'))) (o\P s (k,p)\0^ 
+ / (exp(ifci(z - z )) - exp(ikx(z + z ))) 

J hi 

x / (exp(ik 2 (z - z')) -exp(ik 2 (z + z'^ulolP^h^) 

x (/ - u J (I - exp(2ikzo)P s (k,p))- 1 P s (k 2 ,p)\Oy (54) 

The procedure of Sec I VII which we used to derive Eq.flHUJ) starting with (}4*%|) guarantees that 
in expressions like J k (exp(ik(z — z )) — exp(ik(z + z ))) P s (k, N) entering (jSljl . where k in 
P s (k, N) takes both positive and negative values, the Dirichlet boundary condition at z = 
is correctly taken into account along the contour length of the polymer. 

The poles of the partition function, which are the zero points of the determinant 

det(J -u 1(1- exp(2ikz )P s {k,p)) = (55) 

•i k 

gives the localization energy of adsorbed polymer. The truncation of P s (k,p) by a finite 
order matrix permits to study the adsorption of a semiflexible polymer in a week surface 
potential. Truncation with matrix of order n gives the eigenvalue condition from (|55|) as 
a polynomial of nth degree in powers of u. Due to the same reasons as for adsorption in 
a symmetric potential we will evaluate (|55J) using truncations with even size matrices. We 
will study here only the effect of polymer stiffness on the threshold value of the strength 
of the adsorbing potential u, i.e. Eq. (j55|) for p = 0. For a flexible polymer the threshold 
value of the potential strength u c , such that for u < u c the polymer is delocalized, obeys 
the condition 1 — 3u c zq = 0. The computation using truncation with 4x4 matrices at 
the value z = 1/3 gives the critical value of the localization strength u c = 0.3285, which is 
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smaller than u c = 1 for flexible polymer, and shows that the semiflexible polymer adsorbs 
easier than the flexible one. 



VIII. SEMIFLEXIBLE POLYMERS WITH MONOMER-MONOMER INTERAC- 
TIONS 

The formalism developed in this work enables one to study semiflexible polymers with 
monomer-monomer interactions which are described by the energy U(r(si) — r(s 2 )) (in units 
of ksT). We will restrict ourselves to one polymer, and similar to corresponding treatment 



for flexible polymers 



will consider the correlation function 

•iV rN 



r(*i))) ) , (56) 



G(k, k'; N) = (exp(-tkr(N) - zk'r(O) ~\ J ds 2 d Sl U(r(s 2 ) 

where the average has to be carried out in accordance with Eq.(|T|). To derive the pertur- 
bation expansion of G(k, k'; N) we expand (|5fi|) in powers of the interaction energy, which 
is supposed to be expanded in Fourier integral, order the integrations over the contour 
length, and express all monomer positions entering (|56j) through the tangents according to 
r(s 2 ) — r ( s i) = j'sl dst(s) and r(N) = r(0) + f Q dst(s). The integration over r(0) gives 
the factor (27r) 3 5( 3 )(k + k'). Similar to Sections IVl and IVlTl we finally arrive at the following 
expression under the average 

/ rN rS2n rsi \ 

( exp(— ik / dst(s)) exp(— iQ n / dst(s))... exp(— iQi / dst(s)) ) , (57) 

\ J Sin J &2n-\ JO I 

where the momenta Q n , Qi are expressed by k and q n , qi using the momentum 
conservation in complete analogy to flexible polymers |4j. Using the representation of (J57|) 
through the eigenstates of the quantum rigid rotator we obtain finally the Laplace transform 
of G(k, k'; N) as 

/ - / ^(g n )...^(g 1 )(0,0|P(k,p)P(Q n ,p)...P(Q 1 ,p)|0,0). (58) 

Jq„ Jq 1 

The use of the relation (|18|) enables one to eliminate the propagators P(Qj,p) in favor of 
square matrices P s (Qi,p). Notice that to obtain from (J58|) the corresponding expression 
for a flexible polymer one should replace all P(q,p) in ()58)1 through l/(q 2 /3 + p), and take 
into account in sums only the term with /j = 0. Thus, the comparison of (J58|) with the cor- 
responding expression for a flexible polymer shows that the perturbation expansion of the 
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correlation function G(k, k'; p) for a semiflexible polymer in powers of the interaction energy 
can be represented by the same graphs as those for flexible polymers. However, the associa- 
tion of the graphs with analytical expressions occurs according to Eq. (J58|) . The perturbation 
expansion given by Eq. (j58|) and its straightforward generalization to many polymers is the 
basis for studies of semiflexible polymers with monomer-monomer interactions. 

IX. CONCLUSION 

To conclude, we have developed the statistical-mechanical theory of semiflexible poly- 
mers based on the connection between the Kratky-Porod model of a semiflexible polymer 
and the quantum rigid rotator in an external homogeneous field, and the treatment of the 
latter using the quantum mechanical propagator method. The examples considered in this 
article show that expressions and relations existing for flexible polymers can be generalized 
in a straightforward way to semiflexible ones. The correspondence is established via the 
replacement of the propagator of the theory of flexible polymer through the matrix in the 
theory of semiflexible polymers 

^l—^il + ikDM^D, (59) 

and consideration of an appropriate matrix element of the matrix expression associated with 
the quantity under consideration. The present method provides also a necessary framework 
to study problems including tangents of polymer configurations. 
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